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Abstract: In the case of the gluon-fusion process in deep-inelastic leptoproduction, 
I explicitly show how to incorporate NLO corrections in a Monte-Carlo event gen- 
erator by a subtraction method. I calculate the parton densities to be used by the 
event generator in terms of MS densities. The method is generalizable. A particular 
motivation for treating the gluon-fusion process is to treat diffractive deep-inelastic 
scattering properly, since in diffractive scattering the gluon density dominates the 
quark densities. I also propose a modified algorithm for treating parton kinematics 
in event generators; the new algorithm results in much simpler formulae for the NLO 
corrections. A disadvantage of the new method is that some of the generated events 
may have negative weights. However, an adjustable cut-off function is present in the 
formalism, and this permits a renormalization-group-like transformation that can be 
used to at least reduce the proportion of events with negative weights. 
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1. Introduction 

In an event generator for leptoproduction^ one has a choice of generating events 
either by using the lowest-order parton model process or using the next-to-leading 
order (NLO) hard scattering matrix elements (for photon-gluon fusion, etc). The 
first case is suitable for the total DIS cross section, where one neglects the NLO 
subprocesses, since they represent order corrections to the basic process. The 
second case is suitable, for example, when two-jet production is to be calculated. 

Ideally, one wants to include both the LO and the NLO terms to get good ac- 
curacy. The problem is particularly acute in diffractive DIS, since the gluon density 

^Of course, exactly similar considerations apply to event generators for other processes. 
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is substantially larger than the quark density. Thus the gluon-induced NLO subpro- 
cess is not necessarily smaller than the LO parton model process. This implies that 
inclusion of the photon-gluon term is mandatory to get a sensible phenomenology. 
However, the initial-state showering associated with the LO process already includes 
part of the photon-gluon process, and it is not at all obvious how the two terms are 
to be combined. 

In this paper, I explain how to consistently use both terms in an event generator 
like LEPTO [0] or RAPGAP that uses the algorithm constructed by Bengtsson 
and Sjostrand for initial-state showering. As a consequence of their method 
for treating the exact parton kinematics, the resulting formulae have a non-linear 
dependence on the parton densities. Therefore, I also propose an alternative leading- 
order algorithm for which the corresponding NLO corrections have a conventional 
structure. 

With the exception of the recent paper by Friberg and Sjostrand previous 
attempts, e.g., @, |[, at incorporating NLO corrections have tended to implement 
them by a reweighting of the events generated by showering from the LO matrix 
elements. In contrast, the subtraction method that I describe involves generating 
two classes of events. One class is made from the LO parton-model process by 
showering the initial and final state quarks, exactly as at present. The second class 
of events is generated by starting with a photon-gluon fusion process, and showering 
the partons, again exactly as at present, but with one exception. The exception is 
that the hard cross section for the photon-gluon fusion subprocess is equipped with a 
subtraction that correctly compensates the double counting between the two classes 
of events; the subtraction removes that part of the photon-gluon fusion term that is 
included in the LO parton model plus showering calculation. My method is rather 
similar to the one outlined by Friberg and Sjostrand in Q]. 

Of course, there is no physical distinction between the two classes of events; 
they populate the same regions of state space, and the events differ only in how 
the program generates them. Indeed the relative contributions of the two classes of 
events can be changed by changing the cut-off on the virtuality of lines in the parton 
showering. A change of this cut-off amounts to a renormalization-group transfor- 
mation, and does not affect the physical cross section, aside from the error due to 
uncalculated higher order corrections. 

There is a contrast with the standard subtraction method used in analytic cal- 
culations. In that method an infinite number of positive weighted events is gen- 
erated from the basic photon-gluon fusion graphs. The subtraction term gives an 
infinite number of events with negative weights and somewhat different kinematics. 
Although the infinities cancel in the integral over the hadronic final states, the in- 
tegrand is totally unsuitable for use in an event generator. The new method does 
the subtractions point-by-point in the integrand, so that a Monte-Carlo integration 
is entirely satisfactory. 
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It is true that a finite number of events witli finite negative weights may be 
generated with my algorithm. This is at least a potential disadvantage, particularly 
if it is desired to compare the results directly with data. I will discuss this problem 
in more detail in Sec. ^ There I will show how to overcome this problem, to a 
large extent, by adjustment of a cut-off function that is present in any Monte-Carlo 
algorithm. I will also explain that if an event generator gives a small number of 
negative weighted events, this is not necessarily a severe disadvantage. 

Since an important immediate application is to deal with the problem of the large 
gluon density in diffractive DIS, I give results for the photon-gluon fusion process. 
This case is technically simple because there are no soft gluon effects and no virtual 
graphs. However, the method is capable of being generalized. 

In Sec. 1^, I summarize the basic algorithm used in the event generator at LO, 
and then compute the corresponding first-order term for the photon-gluon fusion 
process. Then in Sec. I explain the observation of Bengtsson and Sjostrand 
that the intended kinematics for the target hadron are not correctly reproduced. 
After reviewing the implementation of their correction to the kinematics, I obtain the 
corrected first-order cross section, from which I compute the subtracted gluon-fusion 
cross section. This forms the primary quantitative result of this paper: it is intended 
to be directly implemented in an event generator. The parton densities are not in 
the MS scheme, so I show how to relate them to the ordinary MS parton densities. 
Unfortunately the resulting formulae are rather complicated; in fact, they are non- 
linear functionals of the parton densities. So in Sec. I present a new algorithm 
for treating the parton kinematics and show that the resulting NLO corrections are 
simpler and more conventional than with the Bengtsson- Sjostrand algorithm. In 
Sec. ^, I discuss the problem that events with negative weights may be generated, 
and show how they may be at least reduced in number. Finally, I summarize the 
directions for future work in Sec. |[ 



2. Basic Monte- Carlo algorithm 
2.1 Algorithm 

I first review the algorithm |[ used in LEPTO or RAPGAP. Only the first part of 
the algorithm will be relevant for our later discussions:^ 

1. Generate values of x and y (and hence Q) from the LO cross section for DIS: 

da 



dx dy 



KF,{x,Q'), (2.1) 



describe the algorithm for the case of fuUy inclusive DIS. The case of diffractive DIS is handled 
by changing the proton to a Pomeron, by replacing the variable x by (3, etc. 
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with 

K = ^^(l-y + y'/2), (2.2) 



and 

F2 = Y.^l^faM^)- (2.3) 

a 

Here, the sum is over all flavors of quarks and antiquarks, and fa is the parton 
probability density. 

2. Generate a virtuality Q\ for the incoming quark a, a longitudinal momentum 
fraction Zi for the first branching, and an azimuthal angle for this branching. 
The distributions arise from the Sudakov form factor 



fc{x/z,Q'^) 




5a(a:,g_,g,) = exp|-y^^ Q^^^?X T^^-^^^^^T 



Here Q^ax is normally set equal to Q^. The Sudakov form factor is the proba- 
bility that the virtuality of the struck quark is less than Q^. 

3. Iterate the branching for all initial-state and final-state partons until no further 
branchings are possible. 

4. Generate 4-vectors for the momenta of all the generated partons. 
2.2 First-order term in Monte-Carlo 

Our aim in this paper is to calculate the NLO contribution to deep-inelastic scattering 
from gluon-fusion graphs like Fig. |l|. The result is to be accurate for the case that 
the incoming gluon, of momentum p^, has virtuality and transverse momentum small 
compared with Q, and that the intermediate quark, of momentum pi, has a virtuality 
of order Q^. To avoid double counting, it is necessary to subtract the contribution 
in this region that is obtained from the showering algorithm applied to the initial- 
and final-state partons of the LO partonic cross section. 

This contribution is obtained by multiplying the lowest order cross section, from 
Eqs. ( |2.1| - |2.3| ), by the the first order term in the expansion of the Sudakov form 
factor ( |2.4| ) in powers of agiQ"^)- This is made differential in the momenta of the 
particles involved, and the gluon-induced term is selected: 

Kj:el^^C{Qi)P{z,)^f,{x/z„ Q% (2.5) 



dx dy dQl dzi dcf) a " 'in'^Qf z\ ' 

Here the splitting kernel is for g quark + antiquark: P{zi) = Pg^qq{zi) = |(1 — 
2zi + 2zl). Note that because we are doing a strict expansion in powers of ^^(Q^), 
the scale argument of the parton density is Q^. The function C{Ql) is a cut-off 
function that gives the maximum value of Q^. To reproduce the upper limit on the 
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Q'"^ integral in the Sudakov form factor, Eq. (p^, one sets C{Q'l) = 9{Q'^ — Q\). 
However, as explained in the introduction, one can envisage changing the cut-off 
function. 

One might worry that a full Sudakov form factor should 
appear in Eq. (p^.S]), to represent the actual physical suppres- 
sion of the cross section at low Q\. In fact, the Sudakov 
form factor should not be used in this formula, since the 
raison d'etre of (|2.5|) is to be a subtraction term for the 
NLO contribution to the cross section. The unsubtracted 
NLO contribution — Eq. (|3.2CI|) below — has the same singular- 
ity and lacks a Sudakov factor. The subtraction will cancel 
the singularity — see Eq. ( p.21| ) — to leave an NLO term that is 
dominantly in the region of large Q\. Thus a strict expansion 
to lowest order in asiQ"^) is appropriate, and a resummation 
of higher-order terms, such as is represented by the Sudakov 
form factor, is not needed. 




Figure 1: 

gluon fusion. 



Photon- 



3. Bengtsson and Sjostrand algorithm 

Bengtsson and Sjostrand give a prescription for the 4-momenta of the partons, 
and we now apply it to our first-order calculation. In this scheme, the 4-vectors for 
the momenta g, pi and of the virtual photon, the intermediate quark and the 
incoming gluon obey the following requirements: 

1. g'^ is the correct value of the photon's momentum. 

2. The proton is to be moving in the —z direction. 

3. pI = -Ql 

4. (pi + g)2 = (p3 - piY =pI = 0. 

5_ ^ Pi ■ {pi + q) _ 
^ P3-{p3 + q)' 

In the 7*(? center-of-mass frame, we therefore have 

_ Q (1 - Ql/Q') 2., \ 



gyi - - QW 



2giv/^i[i-(i + ^i)g?/g2] 2ziQj/Q' 

-riT, i 



Q i-gf/g^ ^' i-Qi/Q' 
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2giv^zi[l-(l + ^i)Q?/Q2] 2z,Ql/Q 



Q l-Q?/g2 ^ ' 1-Ql/Q 

^^^Qj^^m^il.Or.-l). (3.4, 



Here n^- is a unit transverse vector in the direction defined by the azimuthal angle 
0. The components are in the order (0, transverse, 2;). 

We now transform the cross section in Eq. (p.5|) into convenient variables for a 
hard scattering describing photon-gluon fusion. 

The scattering angle obeys 

The fractional momentum of the gluon is defined in [Q] to be 

p-q zi\ Q^J 

where is the proton's momentum. Up to a correction of order m'^/Q'^, this defini- 
tion agrees with the definition in terms of light-front variables: X3 = (^3 — p^) / (p" — 
p^), since the gluon is on-shell with zero transverse momentum. 

We will also use the inverse transformation, to give Qi and Zi in terms of X3 and 

9: 

Ql /I n\ X3 

- = (l-cos^)-, 

z, = ^-hl-cos9). (3.7) 
X3 2 



The Jacobian of the transformation is 

9(0:3, cos ^) 2x 



Then the cross section is: 

J (BSl) 



dx dy dxs dcos 9 d(j) a " 47?^ 



KT^^l^^Cm (3.9) 



^ -p{^-].{l-cos9)]-fg{x/z,,Q' 



1 — cos 6' \X3 2 J x-i' 

The superscript on the cross section is 'BSl' rather than 'BS'. This anticipates that 
we will be forced to change the algorithm; it is the modified cross section that will 
be denoted by a superscript 'BS'. 
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3.1 Inconsistency in kinematics 



Observe that in Eq. (3^) the fractional momentum argument of the gluon density is 
xj Z\ rather than the actual fractional momentum of the gluon, X3. 

This is a symptom of the inconsistency explained by Bengtsson and Sjostrand [Q]. 
The problem can be seen quite dramatically by applying our first-order calculation 
to the case of a gluon target. In that case the density of gluons is a delta-function: 
/g(x3, QX) = 6{x3 — 1), and the generated value of zi is x, from Eq. ([^.51). But then 
the definitions of the parton 4-momenta imply that the fractional momentum of the 
gluon is 0:3 = 1 — Qi/Q'^ — see Eq. ( |3.6|) — instead of the correct value of X3 = 1. 

3.2 Parton model 

What has gone wrong can be explained by examining the derivation of the parton- 
model formula. We consider graphs for deep inelastic scattering that have the form of 
Fig. 1^. For our purposes we can ignore polarization effects, so that the contribution 
of this graph to the structure function is 



F = — 
2n 



(2vr) 



-L{p,p,)Hiq',pipi,)Uipi,] 



(3.10) 



Here L represents the lower part of the graph, H represents the hard scattering, and 
U represents the upper part of the graph. The overall factor in the definition of 
the structure function ensures that it obeys Bjorken scaling, and the overall factor 
1/2tt gives it the standard normalization 




Figure 2: The parton model for deep inelastic scattering is derived from cut graphs of 
this form. 

We now use light-front coordinates [defined by = {V^ ,V~ ,v±), with = 



= (mV2p", p-, 0±) , 
= [Qy2xp~, -xp^, Ox 
Pi = (Pi, ^P~, Pi 



Pi' = \ {Q'^ + xm^)/2xp , {i-x)p , p 



1,-L 



(3.11) 
(3.12) 
(3.13) 
(3.14) 
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The parton-model approximation to the graph is obtained by making approxi- 
mations that are appropriate if the incoming and outgoing quarks, pi and p[ have 
small transverse momenta and virtualities, relative to Q. To make this quantitative, 
we let the magnitudes of plj_, p\, and p\i be M^. Then ^ = x + OiM"^ jQ^). Up to 
a power-suppressed correction, we can replace H by its value with massless quarks, 
if (g^, 0,0), and we can replace the value of p\ = ^p^ in L by xp^ . A shift in the 
integration over ^ then gives the factorized form: 



xp 



dpi (i^pi,_L 



{pt,xp ,Pi,±)) 



/(x)ii(g2,0,0)D + O(MVg' 



i7(g2,0,0) 



dp\, 



dU{pl] 



(3.15) 



Here the parton density /(x) is defined by the usual light-front operator. The frag- 
mentation function D is the integral over the discontinuity of a propagator, so that 
it is equal to unity if the integral is convergent. 

The above derivation is exactly correct in a super-renormalizable field theory, for 
then all the integrals over virtualities are convergent in the approximated integral. 
In a real theory, like QCD, there are ultra-violet divergences that need to be renor- 
malized. The correct physics is obtained in the factorization theorem which shows 
that the above result is valid up to higher order corrections in ^^(Q^)- 

3.3 Correction to the kinematics of initial-state showering 



In the parton-model formula ( p.l5D , the hard scattering has been replaced by a 
massless approximation with zero transverse momentum. Moreover, in the usual 
way of estimating the graph, the quark entering the fragmentation function is given 
zero transverse momentum, and it is given a minus component of momentum equal 
to — x)p~ . For a calculation of the inclusive cross section this use of approximated 
kinematics is correct, since the hadronic final-state is integrated over. But for the 
exclusive calculation, as in an event generator, the use of approximated kinematics 
is wrong, and it results in the inconsistent kinematics noted in Sec. 

The inconsistency is in the definitions of the variable specifying the longitudinal 
momentum of the quark pi. Bengtsson and Sjostrand effectively^ define this 
variable by 

Pi ■ {Pi + q) 



Xi 



p-q 



(3.16) 



■'They actually define a variable z for the relative longitudinal momenta of neighboring partons 
in the showering. Their definition implies that the longitudinal momentum fraction of the struck 



quark is given by Eq. ( p.l6[) . Note that the formula is not in its most satisfactory form. It would 
appear more consistent to replace the denominator hy p ■ {p + q). This results in small changes, of 
relative size xm^/Q^. To keep agreement with the formulae used in I do not make this further 
change. 
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In contrast, the calculation of the Sudakov form factor in Eq. ( |2.4| ) assumes that 
the fractional momentum of the struck quark equals the Bjorken variable x, an 
assumption that is only valid in the limit that transverse momenta, masses and 
virtualities are negligible with respect to Q. 

Hence the inconsistency arises because we have specified the fractional momen- 
tum of pi in two ways. On the one hand the momentum is determined once one 
has specified the virtualities of all the lines, the zi variable for the first branching, 
and the fractional momentum of the incoming parton. On the other hand we asserted 
that the fractional momentum of pi is the parton model value, i.e., that it equals the 
Bjorken variable x. There are several ways one can correct the inconsistency: 

• Change the value of the fractional momentum of pi to the correct value. This is 
the prescription of Bengtsson and Sjostrand 0. We will use it for the remainder 
of this section. It is implemented in LEPTO and RAPGAP, and it requires a 
modified Sudakov form factor. 

• Only generate the value of Bjorken x after the values of the fractional momenta 
are generated. In this approach one cannot choose to generate events with a 
specified value of x. 

• Remove the requirement on the fractional momentum of pi , while keeping the 
unmodified Sudakov form factor. Here one generates the fractional momentum 
variables for all lines except pi, and then determines pi without putting any 
explicit requirement on its fractional momentum. This is the prescription I will 
implement in Sec. |[ 

I now explain the Bengtsson and Sjostrand method. They first make an explicit 
definition equivalent to Eq. (|3.16|) : 



xi=x 1+'-^^^—^]. (3.17) 




Here, m^, = p^, and Qf = —p\. They then write a modified Sudakov form factor: 

(3.18) 

At first sight, this modified Sudakov factor is not satisfactory, since it depends on 
the virtuality of the outgoing jet. This results in a non-factorizing cross section, 
whereas one usually assumes that the showering on different lines can be performed 
independently. In fact, there is no problem as regards the calculation. One simply 
has to perform the generation of the virtuality of the outgoing quark before one 
performs the initial-state showering. A practical alternative is to use the original 
Sudakov factor to generate events and then to reweight the events by the ratio of the 
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form factors. All the necessary variables to define the modified Sudakov factor ( p.l^ ) 
are available when the reweighting is done. An actual reweighting may be done, or 
the generated events may be vetoed with the appropriate probability. 

As we will see, the method also results in formulae for NLO cross sections that 
are non-linear in the parton densities. Similarly, the relations between parton den- 
sities for the event generator and in the MS scheme are non-linear. Although the 
complication in the formulae must be regarded as a disadvantage, the disadvantage 
is not fatal. 

A more fundamental concern is that the dependence of the Sudakov form factor 
on the virtuality of the final-state parton may signal that the form factor is not uni- 
versal between different processes. Since the dependence is only on a single variable, 
one might hope that the problem is tractable. 

3.4 Modification to first-order gluon-fusion 

The modification in the Sudakov factor entails two changes in the cross section 
formula ( |3.9| ). One is that f{x/zi) must be replaced by /(xs), and the other is 
that a factor fa{x)/fa{xi) must be inserted, where Xi = x{l — Q\/Q'^)- Thus the 
first-order cross section implemented in the event generator is 



dx dy dx^dcosO dcj) ^ ^ Att"^ xs fai^i] 



X 



X pf^_i(i-cos^)V (3.19) 

l-cos^ VX3 2^ JJ^ K J 

where the fractional momentum of the gluon in the gluon density is now the same 
as in the calculation of its 4-momentum. 

The extra factor of fa{x)/fa{xi) looks quite unusual. Nevertheless it does repre- 
sent what is actually implemented^ in the code for LEPTO and RAPGAP, following 
the prescription of 0]. This factor does go to unity when Qi ^ Q^, as is necessary if 
the showering is to be correct in the collinear limit. However our purpose is to derive 
the NLO correction to the event generator, and for that we must use the actually 
implemented cross section when is of order Q^. 

3.5 Photon-gluon fusion with subtraction 

From standard references (e.g., 0), we find that the unsubtracted photon-gluon 
fusion contribution to F2 gives a cross section 

C?0"unsubtracted(-^2 part) ^ 2'^s{Q'^)x 2\ nn\ 

^ K 2^ .2 — /9(^3,Q ) X (3.20) 



dx dy dx'i dcos 9 dch ^ a ^^2 

'J T quarks a 

X [piz 



quarks a 

1 



1 — COS 9 1 + cos 9 



^Private communication from H. Jung and T. Sjostrand. 
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where z = x/x^. Notice that the sum is over quark, but not antiquark flavors. 

We must now subtract the 0{as) term associated with the showering, Eq. ( |3.19| ), 
for both quarks and antiquarks. The antiquark term can be treated as a quark term 
with 9 replaced by tt ~ 9. This gives 

^4ard(^2 part) oOsiO^ X ^ 



dx dv dx-i dcos 9 dS z-j a rj. 



quarks a 
1 



X 



+ 



- COS 9 
1 



1 + cos 6^ 



P{z) - CH) ^^p(.-l(l-cos6 

fa{Xlt) V 2 

P{z) - C{-u) ^^P(z -hi + COSi 

fa{Xiu) V 2 




(3.21) 

Here again z = x/x^, and P{z) = |(1 — 2z + 2z'^). The virtualities in the cut-off 
function are —t = Q^(l — cos^^)x3/2x and —u = Q^{1 + cos9)xz/2x] these are used 
to give xu = x{l + t/Q"^) and xiu = x{l + u/Q^). 

Observe how the redefinition of the kinematics of the parton showers has resulted 
in a formula for the NLO correction that is non-linear in the parton densities. In 
principle this is not incorrect. However, the structure of the formula is much different 
to what one usual deals with, and is more complicated. 

For completeness, here follows the corresponding contribution to the cross section 
that results from the photon-gluon fusion part of Fl: 

da\,^rd{FL part) _ Ky"^ ^ ^^ ^.(Q^) x ^ g^) 2^(1 z) (3 22) 
dxdydx3dcos9d(j) 2 - 2y + ^^^^ ^ " Att'^ ^ ^' 



The complete gluon-fusion contribution to the cross section is the sum of (|3.21|) and 

(PD. 



In the above calculations, there is an arbitrary cut-off function C{Qi). To im- 
plement an event generator with the standard Sudakov form factor, one should set 
C{Ql) = 9{Q'^ — Ql). However, the derivation of the leading order formalism does 
not require this choice; the purpose of the derivation is only to obtain the correct 
cross section when Q\ <^ . Another choice of the cut-off function which is unity 
at small Q\ would also be valid. 

Now the aim of computing higher-order corrections, such as Eq. (|3.21|) is to 
obtain the correct cross section when large virtualities are involved. So any change in 
the cut-off function is compensated by the corresponding changes in the subtraction 
terms, up to errors of yet higher order in q;s(Q^). So there is a kind of renormaliz at ion- 
group invariance under changes in the cut-off function. If we were able to calculate 
the cross section to all orders in a^, then the cross section would be exactly invariant 
under changes in C{Q\). 
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Although the cut-off function is arbitrary, it is not completely arbitrary if we are 
to do useful calculations. We must choose it so that higher-order corrections, such as 
Eq. (|3.21|) , are not excessively large. Thus the standard case C{QD = 0{Q^ — Qf) is 
a simple rational choice. However, if one examines the ranges of virtualities actually 
involved, one may well find that some other choice is better. 

Indeed, there is no need to restrict one's attention to sharp cut-offs. A smoother 
function, like 



might even be better, since then one could for example choose it to track the typical 
behavior of a higher-order correction. In this formula A is an adjustable parameter. // 
the cut- ojf function is changed, then one must also redefine the Sudakov form factor. 
Such a change is likely to be particularly useful in connection with the problem that 
Eq. ( p.21| ) does not automatically give a positive cross section. 

This will be discussed further in Sec. |^, where I will show that the use of a 
suitable cut-off function can at least substantially reduce the number of negative- 
weighted events, and where I will also discuss the extent how the use of a formula 
that contains a negative cross section is not necessarily incompatible with its use in 
an event generator. 

3.6 Comparison with MS scheme 

Once one has a systematic algorithm for treating NLO corrections in event generators, 
it is both possible and necessary to answer the question of what scheme is being used 
for the parton densities. The numerical values of the parton densities can be changed 
at order as and beyond by a change of scheme. Thus it is necessary to know the 
actual scheme being used in order that the cross section is actually known to the 
claimed accuracy of a^. Moreover, standard fits to parton densities are typically 
made in the MS scheme; it is necessary to translate these to the scheme used in the 
event generator. 

It does not seem to me to be possible to modify the Monte-Carlo algorithm 
to use MS densities directly, since the algorithm exphcitly uses the dependence on 
parton 4-momentum of the parton correlation function in the target. In contrast, the 
definition of the MS densities involves an integral over all values of the transverse 
momentum and virtuality, followed by a renormalization of the consequent ultra- 
violet divergence. So one must resign oneself to the fact that the event generator 
uses parton densities that are in effect tailored to its algorithm. 

A fundamental way of approaching this issue would be to deduce, from derivation 
of the Monte-Carlo algorithm, the operators that define its parton densities. 

But here I will take a lower-level approach, which is to integrate the Monte-Carlo 
cross section over the hadronic final states and then to require that the result be the 




(3.23) 
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same as the cross section in the standard factorization approach with MS parton 
densities. 

In this approach one must be concerned that the relation between the parton 
densities might be process dependent. In fact, the process dependence cancels. For 
example, the same unsubtracted photon-gluon cross-section Eq. ( p.20| ) occurs in both 
the Monte-Carlo calculation and the MS calculation, and so it cancels out when one 
subtracts the two formulae for the same physical cross section to obtain the relation 
between the parton densities. 

First we take the formula for the structure function with MS parton densities [0]. 
Now our aim is to compute a process-independent relation between parton densities, 
and we want a result for each flavor of quark and not just the combination that 
appears in the usual electromagnetic F2. So it is convenient to replace the actual 
cross section by one in which the photon couples only to one flavor of quark, with 
unit coupling: 



F,"(x,Q^)=x/r)(x,/i2^ 



zvr Jx X3 



Q\l-z) 1 



P{z) In - - + 4z(l - z) 

fj,^z 2 

+ first-order quark terms -|- 0{a1), (3.24) 



where z = x/x^, the same as in the formula for the hard scattering. This must equal 
the same structure function given by the Monte-Carlo calculation. From the order 
as calculation in Eq. ( p.21[ ), we find 



2-n Jx " J -I X3 



+ C dx^ t dcose—fJxs,Q^) X 

Jx J-1 X3 



1 



COS 6 



-^'-"iSK^-5(^-"- 



+ first-order quark terms + O(a^). (3.25) 

Here —t = Q'^{l—cos9)xs/2x, andxi = x— |x3(l— cos^^). Because of the dependence 
on fa{xi), it is not possible to do the integral over cos 6' analytically. The superscript 
on the quark density indicates that it is in the scheme appropriate for the Bengtsson- 
Sjostrand algorithm. 

It follows immediately that 

xe^\x,Q')=xfi^'\x,f,') + 



, I rfa:3-/r)(x3,/i^) 

2tt Jx X3 ^ 
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xlP(z) In 



+ z{l-z) 



1 dcos 9 



P(z) 



-1 1 — cos^^ 
+ first-order quark terms + 0(0;^). 



(3.26) 



This is clearly a rather unpleasant formula. However it is necessitated by the algo- 
rithm currently used in the event generators. 



4. New algorithm for parton kinematics 



The formulae generated in the preceding sections are rather complicated and have 
non-linear dependence on the parton densities. This is a result of the particular 
choice made in [Q for handling the kinematics of the off-shell struck quark. So I now 
propose an alternative algorithm which will result in much more pleasant properties 
for the cross section and the parton densities. 

Consider first how the 4-momenta of the partons are gen- 
erated. In Fig. 1^ is shown a general shower. The showering 
algorithm generates a specification of the partons. The spec- 
ification consists of (a) the flavor of each parton, (b) the vir- 
tuality of each parton, (c) the fractional momentum Z2n+i for 
each branching, and (d) an azimuthal angle for the transverse 
momentum of each branching. The fractional momentum of 
each space-like line is specified as 



X 



=0 



(4.1) 



so that X3 = x/zi, X5 = x/{ziZ3), etc. In reconstructing the 
4-momenta, Bengtsson and Sjostrand tell us to use 



Pi 



P3 



P2n+1 ■ {P2n+1 + Q) 



p-q 



(4.2) 



Figure 3: Generic 
initial-state shower. 



If we apply these equations for all n from upwards, we get the previously explained 
inconsistency. The resolution proposed by Bengtsson and Sjostrand was to change 
the formula, ( ^4.1| ); they made a consistent definition of a fractional momentum X2n+i 

QD/Q"^], and they used the 



[m 



/2 



as the right-hand side of Eq. ( |4.1| ) times [1 -|- 
modified Sudakov form factor ( p.l8| ). 

An alternative which generates simpler formulae is as follows: 



1. Generate virtualities for the partons and fractional momenta for the branchings 
from the unmodified Sudakov form factor Eq. ( |2.4| ). 
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2. Use Eq. ( [4.1|) to compute the fractional momenta, except in the case n = 0, 
i.e., except for computing xi. 



3. Similarly, use Eq. ([4. 21) for all P2n+i, except for pi. 



4. The momentum pi is now completely determined, given the virtualities and 
the azimuthal angles. 

The inconsistency is eliminated because we have dropped one equation, for the lon- 
gitudinal momentum of pi. Bengtsson and Sjostrand, in contrast, changed the com- 
putation of every X2n+i so that the two conditions for pi are consistent. 

The price of the new procedure is that the definition of the longitudinal momen- 
tum fraction, Eq. (O) is not applied universally. The reward, as we will see, is that 



the formulae for the cross sections etc. are now of a conventional form that is linear 
in the parton densities and that is therefore considerably simpler. 

Both procedures are correct; they simply represent alternative prescriptions for 
computing the 4-momenta of the partons from the scalar values generated by the 
algorithm. 

4.1 First-order term for showering 

Since the showers are generated from the unmodified Sudakov factor, Eq. (|2.4|) , the 
cross section for the showering at order as is given by the unmodified Eq. (|2.5|) . The 
parton momenta are given by simpler formulae than Eqs. (|3.1H3.4|):^ 



= I ^ (1 - 2z,, Ot, 1) , (4.3) 



QJl - zi 

pt, = ^ ^ (1, nrsing, cos^), (4.4) 
2\M 



QJl - zi 

p^ = ^ ^ (1, -riTsin^, -cos^), (4.5) 



2^zi 

Q_ 

2^z^{l-zi) 



= , ^ , (1, Ot, -1) . (4.6) 



Here, we now have 



so that 



sin 9 



cose = l-^, (4.7) 
2Qi ' " 



^In the list of requirements above Eq. ( |3.lD , we only need to modify the formula for zi, which 
now becomes zi = Q^/2p3 • [p^ + q). 
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The change of variables is now to cos^^ defined by Eq. ([4.7|) and to X3 = xjzx. 
This gives 



^ (now) 
shower 



dxdydxsdcosOdcp ^ ^tt^ 1 — cos^^ kx^J X3 

4.2 NLO term with subtraction 

The hard scattering cross section with its subtraction is changed to 

'h^\^2 part) ^ y ^2 O'siQ'^ 
dx dy dx'idcosO dd) , " An^ 



C^f^hard (^2 part) ^ 2"£(Q!)a:.. f . . r.^ 

~ K 1^ .0 — /g(^3,Q) (4.10) 



quarks a 

[ P(^)[l-C(-t)] P(.)[l-C(-»)] 1 1 
\ 1-COS0 + l + cos^ 2+^"^^^ ^^/' 

where ^ = x/x^, —t = Q^{1 — cos6')x3/2x and —u = Q^{1 + cos6')x3/2x. Notice the 
considerable simplification compared with the previous scheme. 



4.3 Relation to MS parton densities 

From Eq. ([4.10|) it follows that the structure function F2 is 

+ first-order quark terms + O(a^), (4.11) 



P{z) \n{l/z) ~^ + 3z{l- z) 



where again z = x/x^. Comparison with the structure function expressed in terms 
of the MS densities, Eq. { ^^ , gives 

x/i-)(x,Q2) = x/p)(x,/i2) 



.1 ^ _ r O'^(l-z) 

PU) In ^ ^\ ' +z{l-z) 



+ ^ f'dxs-C'\x„^,^) 
2n Jx X3 ^ 



-|-first-order quark terms + 0{a1). (4.12) 
This is clearly much more tractable than the previous version. 

5. Negative weighted events 



The subtracted NLO cross sections Eqs. (|3.21| ) and (|4.1CI| ) are not automatically 



positive. This is not unphysical, since these formulae only give a component of the 
cross-section. However it is a problem for an event generator when one generates 
separate classes of LO and NLO events. 

The problem is not necessarily fatal. For example, suppose one wishes to make 
a histogram of a differential cross section. Then one can choose to generate weighted 
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events (which if necessary can be passed through a detector simulation, etc). The 
contents of each bin of the histogram are obtained as a weighted sum of events, and 
there is no need to require that individual events have positive weights. 

What would cause a real disaster would be to have a relatively small final answer 
being obtained by a cancellation of large numbers of positive-weighted events and 
large numbers of negative-weighted events. This is not the case for the algorithm 
proposed in this paper, but it is the case for some kinds of analytic calculation. 

Nevertheless, one often prefers to be able to generate unweighted events, for then 
one can genuinely simulate real experimental data. This can be done on the basis of 
generation of weighted events, by standard Monte-Carlo techniques, but only if all 
the weights are positive. 

Luckily there is freedom in the proposed method to reduce the number of negative- 
weighted events, by the choice of the cut-off function C{Ql). Not only can the po- 
sition of the cut-off be changed from its standard value Q^, but the shape of the 
function can be changed, for example to the form in Eq. ( p.23| ). Appropriate changes 
in the showering algorithm will be needed, as explained in Sec. p.5| . At low Ql, 
the cut-off function goes to unity, so that collinear limits are unchanged. But the 
smoother form of cut-off reduces the amount that is subtracted from the bare cross 
section at large t or m in Eqs. ( p.21| ) and ( |4.10| ). For a suitable cut-off function, the 



result is therefore to reduce the number of negative-weighted events, if not eliminate 
them completely. 

Even if one cannot eliminate absolutely all the negative-weighted events, it should 
be enough for practical purposes if only a small fraction of the generated events have 
negative weights. 

6. Future work 

In the context of diffractive deep-inelastic scattering, the NLO corrections calculated 
in this paper are clearly the most urgently needed. Because the gluon density is 
substantially larger than the quark densities, these particular corrections are not 
numerically suppressed in the way that would otherwise be expected of NLO correc- 
tions. The first new result is the calculation of the gluon-induced contribution to the 
hard scattering [Eqs. ( p.21[ ) and ( ^4.101 )]. The second new result is the gluon-induced 



part of the relation between the parton densities to be used in the event generator 
and the MS parton densities [Eqs. and (|4l2|) ]. 



We saw that within the conventional scheme for performing initial-state show- 
ering, the NLO corrections are unusually complicated non-linear functionals of the 
parton densities. This is an inevitable consequence of the algorithm chosen in [^ for 
imposing consistent parton kinematics. Therefore, in Sec. ^, I proposed an alterna- 
tive algorithm for computing the parton 4-momenta. The change in algorithm entails 
a change in the scheme that defines the parton densities used in the event generator. 
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Different numerical values of the parton densities are to be used in each scheme, and 
these numerical values differ from the MS parton densities that are obtained in the 
standard global analyses. 

Further calculations and theoretical developments that are needed include: 

1. Calculation at order of the quark-induced subprocesses. 

2. Numerical calculations of the parton densities from the standard fits, which 
use the MS scheme. 

3. Extension of the cross section calculations to include weak boson exchange. 

4. Extension of the methods to handle other processes in hadron-hadron collisions 
and e^e~ annihilation. 

5. Non- leading corrections to the showering. 

6. Investigation of the non-positivity of the NLO term in the cross sections: Can 
it be eliminated by a suitable choice of the cut-off function? How acceptable is 
it if an NLO term is negative? 

7. Handling more general processes will require a proper treatment of the soft 
region. Technically, I see this as the most difficult problem. 

8. The NLO corrections are valid when the virtuality of the intermediate quark, 
Qi, is of order Q^. There is no singularity at Qi = 0, so the contribution from 
Qi ^ is suppressed by a power of Q^. Even so, events are generated for all 
values of Q\, down to Qi = 0. Thus one may need to modify the Monte-Carlo 
algorithm to handle this kinematic region. 

9. One of the treatments of the kinematics of off-shell partons entails a redefinition 
of the Sudakov form factor in a way that appears non-universal. It needs to 
be understood whether this non-universality is genuine or whether it is only 
apparent. 

10. It would be useful to understand the relation of the Monte-Carlo formalism 
to definitions of the parton densities that involve explicit parton transverse 
momentum. Compare Mrenna's work [0. 

11. Would it be useful to find a way of computing the hard scattering coefficients 
with off-shell matrix elements? There are obviously non-trivial issues of gauge 
invariance that would arise. 

12. Which of the algorithms for calculating the parton 4-momenta is better? Is 
there a yet better algorithm? 
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